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I.  INTRODUCTION 

For  the  past  several  years,  planar  explosive  charge  arrays  have  been 
detonated  on  the  earth's  surface  for  the  purpose  of  subjecting  a  test-bed  to 
compressional  stress  while  retaining  uniaxial  strain.  These  High  Explosive 
Simulation  Technique  (HEST)  tests  are  useful  in  providing  dynamic  in  situ 
compressional  data  for  determining  uniaxial  strain  versus  stress  equations  of 
state  for  various  earth  materials.  (Bratton  and  Higgins,  1979  (Ref.  1): 
Swartz,  et  al ,  1981  (Ref.  2)  and  Jackson  and  Zelasko,  1982  (Ref.  3))*  In 
addition,  planar  explosive  arrays  also  have  been  used  to  approximate  the 
overhead  airblast  impulses  produced  from  surface  point  charges.  The  limited 
practical  size  of  a  HEST  test-bed  (Fig.  1)  causes  nonu*.iaxial  strain  effects 
to  be  generated  by  the  finite  test-bed  edges.  Generally,  the  test-bed  edges 
cause  volumetric  stress  elease  waves  (named  P  release  waves)  and  rotational 
or  deviatoric  stress  waves  (S  release  waves)  inside  the  test-bed  (Rinehart  and 
Lucey,  1981  (Ref.  4)).  The  two  release  waves  assist  in  releasing  the  uniaxial 
soil  stress  to  a  natural  state  of  lithostatic  stress.  Beyond  the  edge  of  the 
HEST  boundary,  the  first  recorded  motion  is  compressional  in  nature  and  ’s 
followed  by  the  deviatoric  or  shear  wave. 

Previous  interest  has  generally  been  focused  on  the  short  duration  (tens 
of  millisecond)  initial  compressional  cycle  of  uniaxial  strain  beneath  the 
HEST.  Attempts  have  been  made  In  understanding  the  edge  release  effects  in 
order  that  the  arrival  time  of  perturbations  from  uniaxial  strain  could  be 
estimated.  The  HEST  would  then  be  sized  to  guarantee  uniaxial  flow  for 
specified  periods  of  time,  within  the  test-bed.  As  late  time  data  became 
available,  it  became  evident  that  In  dry  alluvium  the  release  waves  produced 
stresses  and  strains  on  the  tame  order  as  the  initial  compressional  phase.  To 
completely  understand  cause  and  effects  of  soils  and  structures  *n  the  test¬ 
bed,  an  understanding  of  the  release  waves  was  necessary.  Of  particular 
Interest  in  this  study  are  the  late  time  effects  of  the  release  waves  both 
beneath  the  HEST  and  beyond  its  finite  boundaries.  Constructive  uses  of  these 
effects  for  material  modeling  have  been  identified. 

Long  durations  (up  to  1  s)  of  acceleration  and  velocity  time  histories 
from  several  HESTS  fired  in  dry  alluvium  (Fig.  2)  show  the  edge  release 
effects.  In  all  cases,  the  final  resulting  release  signals  lead  to  wave 
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signatures  that  contain  -1  g  dwells  and  rejoin.  The  -1  g  dwells  indicate  that 
the  soil  is  in  free  flight  under  tensional  failure;  the  material  has  spalled 
from  the  original  surface.  Similar  waveforms,  identified  as  spall,  have  been 
observed  in  data  obtained  from  surface  high  explosive  bursts  (Stump  and 
Reinke,  1983  (Ref.  5);  Merkle,  1980  (Ref.  6}  and  Parret,  1978  (Ref.  7)). 

The  waveforms  include  'nitial  compression  followed  by  release  of  the 
stress.  The  data  show  that  stress  release  caused  by  the  edge  of  the  explosive 
cavity  is  so  severe  that  the  soil  frequently  attains  upward  velocity.  For  all 
records  obtained  for  dry  alluvium,  the  upward  velocity  leads  to  material 
failure  and  free  fall  (-1  g  in  acceleration)  ending  in  rejoin  of  the  material. 
This  spall  initiates  in  such  a  manner  that  the  classical  spall  mechanism  of 
compressional  wave  reflection  (Rinehart,  1959  (Ref.  8);  Eisler  and  Chilton, 
1964  (Ref.  9);  Viecelli,  1973  (Ref.  10))  initiating  the  spall,  is  not  sup¬ 
ported  by  the  observations.  Data  and  calculations  suggest  that  large  magni¬ 
tude  shearing  actions  caused  oy  the  HEST  edges  is  the  probable  cause  of  the 
observed  spall. 

The  primary  objective  of  this  work  is  to  identify  and  quantify  the  mecha¬ 
nism  of  spall  initiation  caused  shearing  motions.  The  mechanism  will  be 
addressed  through  wave  ray  tracing  techniques  and  comparison  of  the  data  to 
two  dimensional  finite  difference  calculations.  In  addition,  a  limited 
parameter  study  was  done  with  the  finite  difference  calculations  to  quantity 
effects  caused  by  changes  in  shearing  properties.  Finally,  some  thoughts  are 
presented  on  extending  these  results  to  spall  phenomena  observed  on  high 
explosive  surface  bursts. 
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II.  SPALL  MECHANISM 

Spall  waveforms  are  a  result  of  tensional  failure  of  in  situ  material 
(Rinehart,  1975  (Ref.  11).  Early  workers  (Rinehart,  1959  (Ref.  8);  Eisler  and 
Chilton,  1964  (Ref.  9)  and  Viecelli,  1973  (Ref.  10))  usually  associated  the 
tensional  failure  with  a  tensional  reflection  from  a  free  surface  of  a 
compressional  wave,  although,  Chilton,  et  al  (1966)  (Ref.  12)  and  Parret 
(1978)  (Ref.  7)  attribute  some  spall  records  observed  for  buried  sources  in 
alluvium  to  excessive  compaction  of  the  alluvium  itself. 

Most  early  observations  of  spall  were  obtained  from  tests  in  hard  rock, 
metals,  or  other  materials  where  tensile  failure  results  in  a  finite  failure 
surface.  The  source  produces  a  high  amplitude  compressive  shock  wave  that 
travels  from  the  source  to  a  free  surface.  Upon  reflection  from  the  free  sur¬ 
face,  a  tensile  wave  is  superposed  onto  the  compressive  wave  until  the  total 
stress  exceeds  the  tensile  limit  of  the  material.  The  upper  spalled  layer 
then  goes  into  free  flight  (-1  g  dwell)  and  is  followed  by  a  sharp  rejoin. 

The  wave  travel  time  curves  for  classical  spall  show  initial  compression 
traveling  from  the  source  to  the  free  surface  and  spall  or  tensional  failure 
traveling  from  ,e  free  surface  to  the  source.  Based  upon  this  physics.  Stump 
and  Reinke  (1983)  (Ref.  5)  ha^e  proposed  criteria  upon  which  spall  can  be 
identified.  j 

!  ■ 

i  ■ 

PRIMARY  CRITERIA  j 

a.  Minus  1  g  (allowable  range  -0.5  g  to  -2.0  g)  vertical  accelera¬ 
tion  dwell.  (Assumes  accelerometers  are  calibrated  to  read  zero  g  in  the 
earth's  gravitational  field.) 

b.  Impulsive  material  rejoin  signals  on  all  components. 

c.  No  acceleration  dwells  on  horizontal  components. 

SECONDARY  CRITERIA 

a.  Dwell  times  of  some  duration. 

b.  Amplitudes  of  rejoin. 
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These  criteria  were  used  throughout  this  work  and  are  adequate  for  spall 
identification  regardless  of  the  initiation  process. 

Using  the  above  criteria  Stump  and  Reinke  (1983)  (Ref.  5)  and  Parret 
(1978)  (Ref.  7)  have  identified  spall  that  does  not  appear,  through  arrival 
time  considerations,  to  be  a  result  of  the  classical  phenomenon.  The  spall 
signals  result  from  both  buried  sources  and  from  surface  point  sources.  As  a 
result  of  the  inadequacy  of  the  classical  spall  mechanism  to  predict  their 
observations.  Stump  and  Reinxe  (1983)  (Ref.  5)  have  proposed  four  additional 
mechanisms  which  provide  a  physically  realizable  source  for  the  required  ten- 
sional  wave  and  failure.  These  include: 

1.  Rayleigh/shear  waves  from  surface  sources--a  wave  equation  effect. 

2.  Spherical  or  cylindrical  divergence  of  the  wave  field— a  geometric 
effect. 


3.  Granular  soil  pore-air  expansion  effect— either  a  source  or  material 
property  effect. 

4.  Material  compaction  and  evacuated  cavity— a  material  property  effect. 


Details  of  these  mechanisms  are  found  in  Stump  and  Reinke  (1983)  (Ref.  5) 
and  Parret  (1978)  (Ref.  7).  A  comparison  of  different  spall  mechanisms  is 
shown  in  Figure  2.  The  top  waveform  clearly  shows  the  features  of  classical 
spall.  The  source  is  an  underground  115  kg  sphere  of  high  explosive  buried 
11.8  m  beneath  the  surface.  The  record  is  taken  from  a  station  between  the 
source  and  the  surface  and  shows  the  arrival  of  the  initial  congressional  wave 
forcing  the  material  upward  followed  by  a  sharp  reversal  and  -1  g  dwell.  The 
reversal  indicates  the  arrival  of  the  reflected  tensile  wave  and  material 


failure.  Note  the  sharp  arrivals  of  both  the  initial  compression  and  tensile 
stresses.  Also  note  the  rapid  reduction  of  the  acceleration  to  -1  g.  The 
second  record  in  Figure  2  is  recorded  from  a  surface  bermed  38  kg  explosive 
sphere.  The  record  is  trken  just  beneath  the  surface  at  a  radial  distance  cf 
a  few  meters  from  the  source.  Initial  motion  is  upward  fallowed  by  a  series 
of  arrivals  that  result  in  a  -1  g  condition  labeled  as  spall  (Stump  and  Reinke, 
1983  (Ref.  5)).  The  source  of  the  tensional  failure  is  believed  by  Stump  and 
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in  initiation  of  the  spall  phase.  The  buried  charge  record  has  a  sharp  rever¬ 
sal,  while  the  surface  burst  has  a  slower,  less  distinct  beginning  of  the 
spall.  These  two  records,  while  both  showing  spall,  have  very  different  pro¬ 
posed  mechanisms.  The  remaining  four  records  are  recorded  beneath  various 
planar  HEST  tests  (Table  1).  For  these  records,  the  initial  phase  is  downward 
compression  followed  by  stress  release  and  a  slow  initiation  of  the  spall. 

The  spall  portion  of  the  records  is  quite  similar  to  the  PHG-79-6  data, 
including  the  slow  initiation  of  the  -1  g  dwell,  rather  than  the  sharp  initia¬ 
tion  seen  in  the  buried  shot  shown.  (As  pointed  out  by  Parret  (1978)  (Ref.  7) 
however,  nonclassical  spall  phenomena  has  been  observed  for  two  buried  explo¬ 
sions.)  The  slow  initiation  of  spall,  in  itself,  suggests  that  a  mechanism 
other  than  classical  spall  is  working  in  the  HEST  test-bed. 

For  surface  point  sources,  the  shear  component  of  the  stress  tensor  is  an 
important  source  of  energy  (Aki  and  Richards,  1980  (Ref.  13))  as  compared  to  a 
buried  source.  Stump  and  Reinke  (1983)  (Ref.  5)  showed  that  the  shear  com¬ 
ponent  created  by  a  point  surface  charge  produces  substantial  tensional 
stress,  which  when  exceeding  any  tensile  failure  limit  of  the  soil  would 
result  in  free  fall.  Present  HEST  data  show  spall  in  the  subsurface  region 
below  the  test-bed  having  similar  characteristics  to  a  point  surface  charge, 
suggesting  a  shear  mechanism.  Spall,  as  will  be  seen,  is  also  present  outside 
of  the  test-bed.  Some  finite  difference  calculations  also  suggest  a  shear 
source. 
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III.  DATA  ANALYSIS 

The  data  used' in  this  reprt  were  obtained  from  motion  measurements  from 
HEST  experiments.  The  explos'v*  array  consists  of  an  explosive  and  bead  foam 
mixture  overlaid  with  a  soil  berm.  Upon  detonation,  the  explosive  and  foam 
mix  and  berm  to  provide  a  peak  gas  pressure,  dependent  upon  the  total  charge 
weight,  within  a  cavity  contained  by  the  berm  on  top  and  soil  underneath  the 
cavity.  The  weight  of  the  berm  controls  the  decay  of  the  peak  pressure  or 
pressure  impulse.  With  a  foam  and  explosive  mixture,  peak  obtainable  pres¬ 
sures  are  much  less  than  the  peak  pressures  produced  by  detonation  of  the 
explosives  alone.  For  material  properties  determination.  HESTs  generally  are 
circular  (DH-1,  DISC-1  and  DISC-2);  although  rectangular  or  square  designs 
have  also  been  used  (SIMCAL  III  and  DAT-3)  (Table  1).  Motions  beneath  the 
test-bed  are  usually  measured  with  high  g  accelerometer  packages  having  uni¬ 
axial,  biaxial,  or  triaxial  capabilities.  The  accelerations  are  then  inte¬ 
grated  to  give  velocities  and  displacements.  In  a  uniaxial  strain  environment 
motion  measurements  alone  are  adequate  for  complete  determination  of  both 
notion  and  stresses.  However,  for  the  relatively  late  times  of  interest  for 
spall  observations  the  strains  are  no  longer  uniaxial.  In  this  case  both 
motions  and  stresses  cannot  be  uniquely  determined  from  motion  gages.  For 
complete  determination  of  the  stress  field,  both  normal  and  shear  stress 
measurements  are  required.  The  state  of  the  art  does  not  presently  allow 
these  measurements;  Indeed,  normal  soil  stress  measurements  are  usually  sus¬ 
pect  in  dry  alluvial  material.  As  a  result,  *oil  stress  is  not  a  parameter 
that  can  be  used  quantitatively  in  this  analysis. 

For  a  recording  station  at  Point  A  (Fig.  1),  beneath  the  explosive  cavity 
the  wave  arrivals  can  be  predicted  using  calculations.  The  initial  motion  is 
vertical,  downward  uniaxial  compression.  Directly  following  the  Initial 
compression  is  the  P  wave  caused  by  the  finite  boundary.  Point  B  (Fig.  1),  of 
the  HEST.  Motion  from  the  P  wave  is  upward  and  outward  within  the  test-bed. 
Outside  of  the  test-bed  the  P  wave  is  compressional  in  nature.  The  next 
arrival  is  the  edge  generated  shear  wave,  also  imparting  an  upward  motion,  and 
Inward  particle  velocity  within  the  test-bed  and  down  and  outward  velocities 
outside  of  the  test-bed.  Point  A  is  also  affected  by  the  release  waves  from 
the  boundary  on  the  opposite  side  of  the  HEST.  For  simplicity  the  opposite 
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edge  release  waves  are  considered  to  be  reflections  from  the  center  line  of 
the  test-bed.  The  position  of  Point  A  and  the  various  wave  speeds  determine 
the  relationship  between  the  arrivals  of  the  various  release  waves.  On  the 
center  line,  both  P  waves  and  both  S  waves  arrive  simultaneously.  In  addition 
to  the  primary  P  and  S  waves,  there  are  P  to  S  conversions  that  appear  as 
diffracted  waves,  with  emergent  arrivals  whose  arrival  times  can  only  be 
estimated. 

\ 

Note  that  beneath  the  test-bed,  all  of  the  release  waveforms  have  motions 
that  imply  stress  release  and  upward  motion.  In  other  words,  under  the  explo¬ 
sive,  velocity  reduction  and  return  to  the  lithostatic  stress  state  is  a 
superposition  of  the  cavity  pressure  decay  plus  the  reversal  due  to  the  finite 
size  of  the  HEST. 

Acceleration  waveforms  recorded  approximately  3  m  out  from  the  center 
line  of  DH-1  for  4.0  m  to  13.0  m  depths  (Fig.  3)  are  seen  to  contain  a  multi¬ 
tude  of  arrivals,  peaks  and  troughs.  Based  on  knowledge  o.  the  wave  speeds 
(Table  2)  and  initial  motions,  an  attempt  was  made  to  pick  out  individual 
arrivals.  As  can  oe  seen,  interpretation  of  individual  waveforms  could  vary. 
What  is  important  is  the  fact  that  the  combination  of  the  edge  effects  make  up 
the  release  portion  of  the  soil  stress.  Spall  initiation,  identified  by  the 
arrival  of  the  -1  g  dwell,  occurs  after  the  reflected  shear  release  wave 
arrives.  Initial  motion  of  the  reflected  shear  wave  is  upward;  however,  the 
acceleration  decreases  slowly  and  becomes  negative.  The  negative  acceleration 
appears  to  be  limited  by  the  lack  of  strength  of  the  material  to  -l  g.  If  a 
more  competent  material  was  being  tested.  It  Is  felt  that  the  acceleration 
could  become  more  negative.  Spall  Initiation  has  a  rapid  apparent  vertical 
velocity.  Implying  a  nearly  horizontally  traveling  wave. 

The  fact  that  the  initiation  of  the  spall  is,  at  any  particular  location, 
not  Instantaneous  with  respect  to  time,  as  compared  to  classical  spall 
(Parret,  1978  (Ref.  7);  Eisler  and  Chilton,  1964  (Ref.  9) ) L  but  rather  a  part 
of  a  wave  train  Is  helpful  In  deciding  on  its  cause.  In  addition,  the  fact 
that  the  spall  appears  not  to  travel  vertically  with  a  proper  wave  velocity 
also  helps  to  rule  against  classical  spall.  \ 

Integrated  accelerations  of  longer  time  durations  (Fiqi  4),  for  the  same 
vertical  cross  section  of  DK-1  show  similar  motion.  The  first  motions  are 
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Initially  down,  followed  by  upward  release,  a  velocity  equal  to  /- 1  g  dt  and 
finally  a  rejoin.  During  the  entire  time  the  relief  effects  are  arriving,  the 
velocity  waveform  indicates  a  smooth  release.  The  initiation  of  the  spall  is 
gradual  and  is  observed  from  the  surface  to  a  depth  in  excess  of  the  radius  of 
the  test-bed.  Extrapolation  of  t  '  data  Indicates  that  spall  reaches  to 
perhaps  li  times  the  test-bed  ra  .  Near  the  surface,  the  spall  phase  is 
followed  immediately  by  a  rejoin.  At  depths  greater  than  6  m,  the  rejoin  is 
more  complicated.  Taking  the  record  at  -IS  m,  (Fig.  4),  the  free  fall  is  ter¬ 
minated,  not  by  a  sharp  rejoin  but  return  to  a  constant  downward  velocity  on 
the  order  of  -1.2  m/s.  A  constant  velocity  of  course,  indicates  zero  acceler¬ 
ation.  This  terminal  velocity  lasts  for  approximately  100  ms  and  is  then  ter¬ 
minated  by  a  rejoin.  The  physical  explanation  of  the  terminal  velocity  is 
unclear  and  does  not  always  appear  in  all  of  the  test-bees.  The  rejoin  appar¬ 
ently  starts  at  the  surface  and  works  its  way  down  at  an  apparent  velocity 
varying  from  3400  m/s  to  207  .i/s.  Unfortunately  the  DH-1  nan  only  vertical 
sensing  gages  beneath  the  test-bed.  On  the  other  hand,  SIMCAL  III  with 
limited  vertical  ranges,  does  have  good  coverage  in  the  horizontal  direction 
both  within  and  outside  of  the  test-bed.  On  the  other  hand,  SIMCAL  III,  with 
vertical  ranges,  does  have  good  coverage  In  the  horizontal  direction.  For 
SIMCAL  III,  information  obtained  from  the  limited  vertical  compares  well  with 
the  DH-1  (Figs.  4  and  5).  Horizontal  motion  from  SIMCAL  III  is  predominantly 
outward.  During  the  spall  phase  horizontal  accelerations  drop  to  zero  as 
required  by  the  criteria;  however,  there  exists  some  small  horizontal  veloc¬ 
ity.  This  Indicates  the  material  that  is  spalled  moves  botn  upward  and  out¬ 
ward.  Beyond  the  test-bed,  which  enas  at  a  range  of  11  m,  spall  is  apparent 
and  is  th?  major  vertical  motion. 


Times  of  arrivals  measured  by  knowledge  of  first  motion  directions  were 
picked  and  compared  to  a  simple  theoretical  model  (Fig.  6).  Because  dry  allu¬ 
vium,  which  at  these  high  pressures  Is  known  to  be  highly  nonelastic  and 
nonlinear,  a  word  about  the  theoretical  model  Is  necessary.  Upon  initial 
loading,  the  stress  wave  will  react  elastically  and  travel  with  a  seismic 
velocity  (Vp).  As  stresses  Increase,  the  stress  stra1n  curve  breaks,  giving  a 
substantially  lower  modulus  (Fig.  7).  The  stress  wave  will  travel  at  a 
loading  wave  velocity  ( Vj_)  of  apphoxlmately  0.4  of  Vp  up  to  a  point  where  the 
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air  filled  voids  are  lost  (on  the  order  of  100  to  200  MPa).  At  any  point 
wnere  unloading  and  subsequent  reloading  occurs,  the  stresses  will  travel  at 
an  unloading  velocity  (VUp)  that  is  near  the  original  seismic  velocity.  For 
this  problem,  a  tri linear  material  model  will  be  assumed  where  the  initial 
uniaxial  strain  loading  will  travel  at  Vp,  the  peak  stress  at  V|_  and  the  volu¬ 
metric  relief  effects  at  Vup  made  equal  to  Vp.  Comparison  with  actual  in  situ 
data  (Fig.  7)  indicates  that  this  is  not  a  bad  assumption.  The  shear  modulus 
is  assumed  to  be  directly  proportional  to  the  constrained  modulus  through  a 
Poisson's  ratio  of  0.33.  Note  that  neither  the  model  nor  the  real  physics 
will  allow  any  releaf  effects  to  travel  faster  than  initial  loading  velocity, 
although  the  speeds  can  cause  stress  release  of  the  initial  compression.  The 
simple  model  that  was  used  was: 

Vp  =  700  m/s 
VL  *  280  m/s 
Vup  *  700  m/s 

Vs  *  330  m/s  (Poisson  ratio  =  0.33) 

which  appears  to  fit  the  data  (Fig.  6)  well,  with  the  exception  of  the  shear 
waves.  The  initial  shear  wave  appears  to  travel  at  330  m/s  beneath  the  explo¬ 
sive  cavity  in  highly  disturbed  material.  The  shear  velocity  increases  to 
380  m/s  outside  of  the  test-bed  in  less  disturbed  material. 

Spall  initiation  travels  in  the  horizontal  direction  at  an  apparent 
velocity  of  560  m/s.  Since  this  apparent  velocity  is  less  than  tNe  volumetric 
reloading  P-wave  velocity,  it  must  be  a  shear  wave  traveling  obliquely  to  the 
horizontal.  An  obliquely  traveling  P  release  wave  would  travel  at  a  velocity 
greater  than  the  true  P  release  wave.  For  SIMCAL  III,  at  this  particular 
depth,  the  -1  g  dwell  is  immediately  followed  by  soil  rejoin;  there  is  no  ter¬ 
minal  velocity  noted  as  in  DH-1. 

Figure  5  shows  that  under  the  explosive  array  the  displacement  is  down¬ 
ward,  followed  by  upward  and  outward  vectors  until  the  spall  commences.  The 
outward  displacements,  however,  are  only  from  5  to  10  percent  of  the  vertical. 
Outside  the  test-bed  at  a  range  of  7  m  from  the  edge  of  the  explosive  array 
(Fig.  8),  initial'  displacement  is  outward  and  down.  The  motion  then  becomes 
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approximately  elliptical  and  retrograde  with  respect  to  the  direction  of  prop¬ 
agation  until  spall  initiation,  at  which  point  the  horizontal  velocity  drops 
to  nearly  zero.  The  displacement  vector  suggests  that  spall  is  directly  pre¬ 
ceded  by  a  Rayleigh  surface  wave,  however,  these  motions  were  observed  in  an 
area  that  still  is  undergoing  nonelastic  strains.  In  other  words,  the 
motions  are  still  probably  a  direct  result  of  the  source,  but  are  clearly  in 
the  transition  zone  between  nonelastic  source  effects  and  true  geologic 
influences.  No  data  exist  to  properly  define  the  transition  zone  for  HEST 
tests. 


A  vertical  profile  >f  displacement  time  histories  was  provided  by  DH-1 
(Fig.  9)  by  doubly  integrating  accelerations.  The  number  to  the  right  of  each 
curve  indicates  the  depth  from  which  the  data  are  taken.  Differences  between 
the  data  traces  at  particular  times  are  proportional  to  the  volumetric  strain 
through  the  spatial  derivative.  Initial  motion  of  the  test-bed  is  downward, 
with  peak  strains  occurring  near  the  surface.  Upon  release,  the  material 
remains  crushed  up,  with  little  strain  rebound  as  expected  with  alluvium, 
bpall  occurs  with  the  material  still  under  initial  compaction,  although  the 
material  is  in  tension.  From  the  surface  to  approximately  5  m  in  depth,  the 
material  regains  some  of  its  initial  volume,  however,  below  6  m  the  entire 
test-bed  remains  compacted  during  the  spall  and  fallback  to  the  position 
obtained  during  initial  compaction. 

In  dry  alluvium,  spall,  implying  tensional  failure,  does  not  imply 
bulking  of  the  material.  The  initial  stresses  are  sufficient  to  cause  signi¬ 
ficant  permanent  and  irrecoverable  compaction.  This  phenomenon  Is  material 
dependent  ard  is  not  expected  to  occur  In  materials  which  do  not  contain  large 
amounts  of  air  filled  voids. 

To  see  if  spall  s  common  among  dry  alluvial  sites,  DISC-1,  DISC-2  and 
OAT-3  were  checked  (Table  1).  In  all  cases,  spall  very  similar  to  what  has 
been  described  existed  (Fig.  2).  Individual  specific  material  properties 
varied  from  site  to  site  (Fig.  7);  however,  they  all  exhibited  the  general 
alluvial  bilinear  material  behavior.  All  of  the  tests  were  approximately  the 
same  size,  so  that  an  Idea  of  the  dependence  of  spall  with  respect  to  test 
size  could  not  be  attempted. 
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IV.  CALCULATIONS 

To  obtain  better  insight  into  the  problem  and  to  be  able  to  identify  the 
causes  and  effects  of  the  spall,  several  calculations  were  performed.  The 
problem  deals  with  material  response  to  shock  conditions  in  a  highly  non¬ 
elastic  and  nonlinear  region  of  material  behavior.  In  addition,  tensile 
failure  and  cracking  are  an  important  aspect  of  the  problem.  As  a  result,  no 
attempt  has  been  made  to  model  the  problem  with  the  usual  elastic  closed  form 
solution  (Murphy  1979  (Ref.  14),  and  Day  et  al ,  1983  (Ref.  15)).  Rather  a 
two-dimensional  finite  difference,  general  computer  code  modeling  energy  and 
mass  conservation  was  chosen  (Trulio,  1966  (Ref.  16)  and  Schuster  1982 
(Ref.  17))  to  investigate  the  data. 

The  problem  is  set  up  with  symmetry  about  the  central  axis  of  the  test¬ 
bed.  For  most  of  the  tests,  except  SIMCAL  III  and  DAT-3,  this  was  the  field 
configuration.  To  calculate  square  or  rectangular  test-beas  accurately  a  two- 
dimensional  code  even  employing  plane  strain  would  be  inappropriate  because  of 
the  four  edge  boundaries.  Recent  work  (Shinn  1983  (Ref.  18)),  however,  has 
indicated  that  to  a  first  order  approximation,  an  axisymmetric  calculation 
represents  square  test-beds  adequately.  In  addition,  by  using  a  simple 
geometry  individual  waves  could  be  studied  from  a  particular  boundary  without 
having  to  contend  with  traveling  waves  from  additional  boundaries. 

A  Lagranglan  formulation  was  chosen  to  prevent  mixing  and  loss  of 
material  near  the  areas  of  large  deformations.  The  grid  was  to  slide  on  the 
axis  of  symmetry,  free  at  the  soil -air  Interface  and  transmitting  at  the 
remaining  sides.  The  grid  was  made  large  enough  to  prevent  any  spurious  grid 
boundary  reflections  from  causing  perturbation  In  the  zone  of  interest  for  the 
duration  of  the  problem.  The  explosive  cavity  was  simulated  by  reasonably 
matching  experimental  pressure  wave  forms  from  SIMCAL  III  with: 

p  x  pQe-<*t  +  pa»x  <  7  m  (HEST  boundary) 

P  -  Pa  ,  X  >  7  m 

where 

P„  -  10  MPa 
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Pa  *  Atmospheric  pressure 

500  for  the  plastic  case 

a  * 

50  for  the  elastic  case 
t  =  Time  (in  seconds) 

The  equation  of  state  of  the  material,  modeled  with  the  AFUL  engineering 
model  (Shuster  1982  (Ref.  17)),  fits  the  uniaxial  strain  versus  stress  data 
of  the  material  with  a  series  of  three  straight  lines.  Plastic  flow  in  shear 
is  allowed  when  the  stress  difference  exceeds  a  Drucker-Prager  failure  cri¬ 
terion  given  by: 

£L^»3.  >  C„  *  8  7 

where 

oi  *  Maximum  Principal  Stress 

<j3  Minimum  Principal  Stress 

C0  *  Cohesion 

6  *  0.7 

F  «£>•/ 3 

.  i 

i 

Plastic  flow  Is  modeled  as  nonassociati ve,  keeping  P  constant  when  calculating 
plastic  strain.  A  simple  cracking  model  with  appropriate  volumetric  adjust¬ 
ments  is  cllowed  when  tension  exceeds  tq,  the  tensile  limit.  Although  this 
model  appears  to  be  quite  simple  it  has  successfully  modeled  many  high 
pressure  explosive  events  and  Is  adequate  for  our  purposes. 

1 

Three  basic  calculations  were  accomplished  with  the  plastic-elastic  model 
shown  jin  Figure  7.  Their  basic  difference  was  the  C0  used  for  limiting  the 
plastic  flow.  In  addition,  an  elastic  calculation  was  done.  The  elastic- 
plastli  soil  model  was  chosen  so  that  the  elastic  P  wave  velocity  is  equal  to 
the  unloading  wave  velocity  of  the  plastic  case.  All  of  the  release  phenome¬ 
non  In  the  plastic  calculations  should  appear  to  travel  at  the  same  velocities 
as  In  the  elastic  case.  Comparisons  can  then  be  made  between  the  two  calcula¬ 
tions  with  respect  to  wave  arrivals,  causes  and  effects. 
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The  calculated  velocity  vectors  at  38  ms  into  the  problem,  for  the 
elastic  case  clearly  shows  the  expected  waves  (Fig.  10a).  The  pressure  puise 
was  made  short  so  that  it  would  not  obscure  the  phenomenon.  Individual  arriv¬ 
als  are  labeled  with  direct  waves  referring  to  those  generated  by  the  HEST 
boundary  located  in  the  problem;  reflected  waves  refer  to  those  generated  by 
the  image  boundary  or  reflected  from  ine  center  of  the  asymmetric  problem. 

The  first  motions  are  those  observed  in  tne  data.  The  direct  P-wave  is 
centered  around  a  P-nodal  plane  extending  directly  beneath  tne  boundary  of  the 
HEST.  First  motions  of  this  P-wave  are  reversed  from  beneath  the  HEST  to  the 

outside.  The  action  of  the  shear  waves  rotating  the  material  upward  beneath 

the  explosive  cavity  is  clearly  seen.  The  arrival  of  the  shear  wave  beneath 

the  test-bed  is  oblique  to  the  surface  as  the  data  suggest. 

Figures  10b,  c,  and  d,  show  velocity  vector  plots  for  the  three  plastic 
cases  at  approximately  the  same  times.  The  pressure  pulse  was  increased  in 
impulse  to  replicate  actual  test  data  better.  The  arrival  of  the  individual 
waves  are  not  as  clear  as  in  the  elastic  case.  The  major  difference  in  the 
velocity  field  oetween  the  three  el astic-plastic  calculations  is  changes  in 
the  plastic  flow  surface  through  changes  in  the  modeled  unconfined  c0.  The 
unconfined  cohesion  essentially  dictates  the  allowed  limit  of  stress  difference 
allowed.  Although  simply  modeled,  it  appears  to  be  an  adequate  representation 
of  the  observed  physics. 

The  calculations  indicate  that,  for  material  having  less  cohesion,  the 
shear  wave  becomes  progressively  less  strong  and  is  absent  from  the  calcula¬ 
tion  with  no  cohesion.  For  the  case  with  zero  cohesion,  the  material  under 
the  explosive  fails  to  rebound  substantially,  and  fails  to  obtain  any  upward 
velocity.  Instead,  plastic  flow  dominates  and  the  test-bed  material  flows 
downward  and  outward  from  the  test-bed.  Observations  showing  this  behavior 
have  been  noted  In. large  HEST  tests  fired  over  saturated  clayey  soils  with 
little  cohesion.  As  the  modeled  cohesion  is  increased,  the  material  is 
affected  more  by  the  shear  wave  and  is  seen  to  obtain  upward  velocities  which 
terminate  in  tensional  failure  and  spall.  Spall  then  appears  and  can  be  asso¬ 
ciated  with  shear  wave  action.  The  presence  or  absence  of  spall  is  directly 
related  to  the  modeled  cohesion  of  the  material. 
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With  the  model  employed,  material  cracking  is  noted  (Fig.  11)  both  beneath 
and  outside  the  explosive  cavity.  Initial  compression  causes  shear  failure 
(in  the  material  with  high  cohesion)  or  plastic  flow  (material  with  ^ero 
coehsion),  allowing  material  to  move  both  downward  and  outward.  The  arrival 
of  the  direct  and  reflected  P-  and  S-waves  causes  the  test-bed  to  be  thrown 
upward  ending  in  tensional  failure  or  spall  as  shown  by  the  horizontal  crack 
indicators.  The  shear  wave  appears  to  be  a  significant  factor  in  causing 
spall,  implying  that  unconfined  cohesion  of  the  soil  is  required  for  spall  to 
be  initiated  by  a  shear  wave.  The  region  of  spall  in  the  calculations  extends 
to  approximately  1|  times  the  radius  of  the  explosive  cavity  agreeing  with 
extrapolation  of  the  data.  As  with  the  data,  failure  is  also  observed  outside 
of  the  test-bed.  The  failure  occurs  rapidly,  as  seen  in  the  data,  requiring 
only  a  few  cycles  of  the  calculations. 

Vertical  velocity  time  histories  for  four  cases  of  different  C0  clearly 
show  the  effects  of  the  shear  wave  interaction  (Fig.  12)  on  the  final  velocity 
time  history.  The  time  histories  shown  are  for  the  centerline  of  the  test¬ 
bed  so  that  the  direct  and  reflected  waves  arrive  simultaneously.  The  time 
histories  again  show  what  has  been  described  before.  Spall  initiation,  iden¬ 
tified  by  a  constant  slope  in  the  velocity  field  equal  to  -1  gt,  does  not 
occur  at  initial  shear  arrival  but  is  delayed  until  the  positive  vertical 
velocity  is  sufficient  to  cause  failure.  The  plastic  case  having  no  cohesion 
does  not  appear  to  permit  shear  wave  propagation,  resulting  in  small  velocity 
reversals. 

The  P-wave  nodal  plane  generated  by  the  finite  boundary  is  clearly  seen 
In  the  calculation.  Near  the  surface,  almost  no  P-release  wave  motion  is 
observed.  As  one  goes  deeper,  the  amplitude  Increases.  However,  if  stress 
time  histories  are  plotted,  the  effect  of  the  P-release  wave  Is  clearly  seen. 

The  normal  stresses  are  substantially  reduced  to  zero. 
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V.  DISCUSSiON  AND  CONCLUSION 

The  importance  of  this  work  is  to  uniquely  identify  a  physical  process 
that  produces  a  spall  waveform,  which  is  not  dependent  upon  the  classical 
reflection  of  a  compression  wave.  As  early  as  1966,  Chilton,  et  al ,  (Ref.  12) 
suggested  alternatives  to  the  classical  theory.  A  great  deal  of  interest  was 
shown  when  spall  identified  with  surface  explosions  could  not  be  related  to 
the  classical  explanation  (Merkle,  1980  (Ref.  6);  Stump  and  Reinke,  1983  (Ref. 
5);  Day.^et  al ,  1983  (Ref.  15)  and  Auld  and  Murphy,  1979  (Ref.  19)).  As  a 
result,  four  additional  mechanisms  have  been  proposed  by  Stump  and  Reinke 
(1983)  (Ref.  5),  in  addition  to  the  classical  mechanism.  This  particular  data 
set  is  of  interest  because  the  final  motion  resulting  from  a  dynamic  uniaxial 
compression  test  is  tensional  failure  and  spall.  Clearly,  the  tensional 
failure  cannot  be  caused  by  a  compressional  wave  traveling  upward  to  the  sur¬ 
face  and  reflecting.  This  required  an  alternate  theory  which  allows  large 
velocity  reversals  through  shear  wave  interaction  upon  release  of  the  soil 
stress. 

With  a  HEST-type  loading  it  is  tempting  to  relate  the  observed  spall  with 
a  surface  point  source.  As  one  gets  away  from  the  test-bed,  the  HEST  approxi¬ 
mates  a  point  source.  Aki  and  Richards  (1980,  p  218)  (Ref.  13)  illustrate 
the  generation  of  shear  waves  from  a  surface  point  source  by  using  a  torus 
with  downward  and  outward  rotating  motion.  In  a  sense,  the  circular  HEST  may 
be  considered  as  just  such  a  geometric  array.  The  calculations  indicate  that 
major  shearing  exists  during  this  Initial  compressional  phase  outside  the 
test-bed  as  indicated  by  theory;  however,  it  is  the  release  wave  from  the 
edges  of  the  test-bed  that  cause  rotation  and  upward  motion  that  results  in 
spall.  For  a  point  source  explosive,  the  edge  effects  are  not  as  clearly 
defined  as  for  the  HEST.  However,  at  any  particular  instant  In  time,  the 
overhead  airblast  from  a  surface  burst  does  have  finite  boundaries  associated 
with  the  airblast.  The  finite  boundaries,  although  moving,  would  cause  an 
initiation  of  a  shear  wave.  Thus  spall  associated  with  airblast  can  be  postu¬ 
lated  to  exist.  In  addition  to  the  airblast,  spall  from  bermed  explosives 
without  airblast  are  also  observed  (Fig.  2).  The  spall  mechanism  observed 
with  the  bermed  charges  Is  in  addition  to  the  spall  due  to  the  airblast  and  is 
probably  due  to  the  same  shearing  mechanism. 
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Comparison  of  data  and  calculations,  indicates  that  the  spatial  extent  cf 

t 

the  spall,  and  thus  its  initial  duration!,  or  wave  length,  are  directly  propor¬ 
tional  to  the  size  of  the  HEST.  Both  thf  data  and  calculations  indicate  that 
the  extent  of  spall  in  a  vertical  direction  is  approximately  li  times  the 
radius  of  the  HEST.  Extending  this  to  a  surface  burst  would  imply  that  the 
vertical  extent  of  spall  is  related  somehow  to  the  peak  overpressure  areal 
history.  Since  the  overpressure  scales  as  the  cube  root  of  the  energy  of  the 
point  source,  the  extent  of  the  spall  on  a  surface  burst  could  also  scale  as 
the  cube  root  of  the  yield.  However,  the  tensional  stress  required  to.cause 
spall  would  appear  to  increase  due  to  increased  lithostatic  pressures. 

Another  important  result  from  this  data  is  the  relationship  between  • 
strength  or  C0  of  the  soiland  the  resulting  upward  vertical  velocity  attained 
just  prior  to  the  initiation  of  spall.  Figure  11  qualitatively  descriDes  the 
limits  of  peak  upward  velocities  obtainable  for  various  C0‘s.  No  upward 
velocity  for  soil  modeled  with  zero  cohesion  to  large  upward  motions  for 
elastic  material  were  calculated.  A  ratio  of  peak  upward  velocity  calculated 
compared  to  the  peak  downward  velocity  was  found  to  be  related  to  the  cohesion 
of  the  soil  model  (Fig.  13).  In  addition,  the  qualitative  appearance  of  the 
spall  initiation  is  also  related  to  the  modeled  cohesion.  The  calculations 
indicate  that  the  cohesion  of  0.5  MPa  produced  a  sharper  spall  initiation  than 
one  of  0.3  MPa.  Ratios  of  the  velocities  can  be  made  using  data  which  indi¬ 
cate  that,  for  soils  tested,  a  cohesion  of  0.3  MPa  to  0.5  MPa  is  required  to 
model  the  spall  signature. 

For  dry  alluvial  soil,  these  cohesions  appear  a  bit  large.  On  the  other 
hand,  a  true  in  situ  estimate  of  cohesion  is  one  property  that  remains  almost 
unmeasurable  with  available  techniques  (A.  E.  Jackson,  personal  communica¬ 
tion).  In  addition,  the  materials  for  which  these  estimates  are  made  have 
been  initially  compacted.  The  displacement  time  histories  (Fig.  9)  also  show 
that  tensional  failure  occurs  while  the  material  is  still  under  uniaxial  - 
strain.  This  also  may  affect  the  final  cohesion  of  the  soil,  before  spall. 
Although  the  relationship  between  cohesion  and  spall  signature  is  now  only  an 
observation.  It  may  prove  to  be  an  important  measure  of  soil  cohesion. 

Spall  and  tensional  failure  Is  observed  in  uniaxial  dynamic  compression 
material  properties  tests,  apparently  a  result  of  shear  waves  generated  from 
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the  finite  boundaries  of  the  experimental  test-bed.  Calculations,  both 
elastic  and  elastic-plastic,  using  simple  soil  models  have  helped  to  confirm 
the  phenomenon.  The  spall  signature  is  much  different  from  that  observed  from 
classical  reflection  of  compressive  stress  waves  impinging  upon  a  free  sur¬ 
face.  The  phenomenon  may  be  applicable  to  surface  point  sources.  Finally, 
calculations  indicate  that  the  spall  signature  and  the  ratio  of  peak  upward 
velocity  versus  peak  downward  velocity  may  be  directly  related  to  the  cohesion 
of  the  dry  alluvium. 
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Figure  1.  Test  configuration  showing  typical  test-bed  and  wave  fronts.  (The 
arrows  indicate  the  direction  of  first  motion.) 
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DH-I 

VERTICAL  CROSS  SECTION 


TIME  (mt) 


VWICiL  ACCILt  RATION* 


NOTE:  The  four  release  waves  and  spall 
Initiation  are  shown.  Records  are  trun¬ 
cated  at  spall  initiation.  See  Figure  4 
for  longer  time  durations. 

Figure  3.  Vertical  accelerations  from  DH-1 
for  a  cross  section  3  m  from  the 
centerline  of  the  HEST. 
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TABLE  2.  PARAMETER^  USED  IN  THE  CALCULATIONS 

I 


Calculation 

Type 

7TT 

(m/s)  j 

VL  , 

(m/s) 

VUD 

(ms) 

Lo 

(MPa) 

Poisson's 

Ratio 

1 

Elastic 

800  | 

800 

800 

N/A 

1/3 

2 

Plastic 

550 

333 

800 

0.0 

1/3 

3 

Plastic 

550  | 

333 

800 

0.3 

1/3 

4 

Plastic 

550 

333 

800 

0.5 

1/3 
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Figure  4.  Velocity  time  histories  from 
DH-1.  (The  location  of  the 
cross  section  Is  3  m  from  the 
centerline.) 
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Figure  7.  Uniaxial  strain  in  situ  material  models  estimated  fpr  the 
various  sites.  [Also  shown  are  the  two  (elastic  and 
elastic-plastic)  models  used  in  the  calculation.]  (After 
Ref.  2.1 


Figure  8.  Displacement  map  of  a  superposition  of  both  horozontal 
and  vertical  displacements.  (This  is  for  a  station  7  m 
from  the  edge  of  the  test-bed  at  a  depth  of  5.4  m  for 
SIMCAl  III.) 
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NOTE:  Plots  are  taken  at  the  times  when  the  reflected  shear  waves  appear  directly  under 
the  test-bed.  The  modeled  cohesion  varies  from  0.0  MPa  to  elastic. 

Figure  10.  Velocity  vector  plots  from  the  calculations. 
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Figure  10.  Continued. 
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Figure  11.  Plot  of  the  grid  and  fracture  pattern  at  three 

times  for  the  case  where  C0\  *  0.5  MPa.  [The  crack 
symbols  (dashes  within  the  gpd)  are  orientated 
normal  to  the  least  principal  stresses.] 
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VERTICAL  VELOCITY 
(NOT  SCALED) 


Figure  12.  Velocity  wave  forms  at  the  symmetry  axis  for  four 
calculations. 
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Figure  13.  Comparison  of  peak  particle  velocity  obtained  at  spall 
normalized  by  the  peak  pressure  compared  to  the  modeled 
cohesion. 
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